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The knowledge of quantum phase flow induced under the Weyl's association rule by the evolution 
f~^ , of Heisenberg operators of canonical coordinates and momenta allows to find the evolution of sym- 
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f^ ■ hols of generic Heisenberg operators. The quantum phase flow curves obey the quantum Hamilton's 

^ ' equations and play the role of characteristics. At any fixed level of accuracy of semiclassical expan- 

j^ ' sion, quantum characteristics can be constructed by solving a coupled system of first-order ordinary 

differential equations for quantum trajectories and generalized Jacobi fields. Classical and quantum 

Cn , constraint systems are discussed. The phase-space analytic geometry based on the star-product 

operation can hardly be visualized. The statement "quantum trajectory belongs to a constraint 

cn 

^ ' submanifold" can be changed e.g. to the opposite by a unitary transformation. Some of relations 

in ■ 

>— P , support partly geometric relations of belonging and intersection. Quantum phase flow satisfies the 

("^ . star-composition law and preserves hamiltonian and constraint star-functions. 
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P3 ■ I. INTRODUCTION 



permits formulation of quantum 



7^ ' The star-product operation introduced by Groenewold for phase-space functions [1] 

mechanics in phase space. It uses the Weyl's association rule 2|, y] to establish one-to-one correspondence between 
Vh ' phase-space functions and operators in the Hilbert space. The Wigner function [J| appears as the Weyl's symbol 



iy, 



governs 



of the density matrix. The skew-symmetric part of the star- product, known as the Moyal bracket 

the evolution of symbols of Heisenberg operators. Refined formulation of the Weyl's association rule is proposed 

by Stratonovich [7L The Weyl's association rule, star-product technique, star- functions, and some applications are 



reviewed in Refs. 



m 



yiiiU- 



A one-parameter group of unitary transformations in the Hilbert space 

ii = exp(--^i3r), (LI) 

with Sj being Hamiltonian, corresponds to a one-parameter group of canonical transformations in the classical theory 
[3|, lla, llll J although canonical transformations provide a broader framework [18J, |l9[ . 

Weyl's symbols of time dependent Heisenberg operators of canonical coordinates and momenta induce quantum 
phase flow. Osborn and Molzahn 20] construct quantum Hamilton's equations which determine quantum phase 



flow and analyze the semiclassical expansion for unconstrained quantum-mechanical systems. An earlier attempt to 
approach these problems is undertaken in Ref. 21|. 



The infinitesimal transformations induced by the evolution operator (jl.ip in phase space coincide with the infinites- 
imal canonical transformations induced by the corresponding Hamiltonian function [3|, lla . Il7| . The quantum and 
classical finite transformations are, however, distinct in general, since the star- and dot-products [56[ as multiplication 
operations of group elements in quantum and classical theories do not coincide. The quantum phase flow curves are 
distinct from the classical phase-space trajectories. This fact is not well understood (see e.g. Refs. 



21 



m 



Osborn and Molzahn [20.] made important observation that quantum trajectories in unconstrained systems can be 
viewed as a "basis" to represent the evolution of quantum observables. 

Such a property is usually assigned to characteristics appearing in a standard technique for solving first-order partial 
differential equations (PDE). The well known example is the classical Liouville equation 

|:/(C,r) = {/(e,r),H(e)}. (1.2) 

This equation is solved in terms of characteristic lines which are solutions of classical Hamilton's equations 

^c\tr)^{CMO}k=cii..r) (1.3) 

with initial conditions c'(^, 0) = ^*. Equations (jl.3[) are characteristic equations. They represent a system of first-order 
ordinary differential equations (ODE) for canonical variables. Physical observables /(^,t) evolve according to 

/(e,r)=/(c(C,r),0). (1.4) 

It is remarkable that despite quantum Liouville equation is an infinite-order PDE its solutions are expressed in 
terms of solutions of the quantum Hamilton's equations which are infinite-order PDE also. 

A technical advantage in using the method of characteristics in quantum mechanics stems from the fact that to any 
fixed order of the semiclassical expansion the quantum Hamilton's equations can be viewed as a coupled system of first- 
order ODE for quantum trajectories and generalized Jacobi fields obeying certain initial conditions. The evolution can 
be considered, respectively, as going along a trajectory in an extended phase space endowed with auxiliary degrees of 
freedom ascribed to generalized Jacobi fields. The evolution problem can be solved e.g. numerically applying efficient 
ODE integrators. 

Quantum characteristics can be useful, in particular, for solving numerically many-body potential scattering prob- 
lems by semiclassical expansion of star-functions around their classical values with subsequent integration over the 
initial-state Wiener function. Among possible applications are transport models in quantum chemistry and heavy-ion 



collisions 



23 



24 



25| where particle trajectories remain striking but an intuitive feature. 

A covariant extensions of quantum molecular dynamics (QMD) transport models [26l . |27| is based on the Poincare 
invariant constrained Hamiltonian dynamics j28| . 

We show, in particular, that quantum trajectories exist and make physical sense in the constraint quantum systems 
also and play an important role similar to that in the quantum unconstrained systems. 



The paper is organized as follows: In Sects. II and III, characteristics of unconstraint classical and quantum systems 
are discussed. Sects. IV and V are devoted to properties of characteristics of constraint classical and quantum systems. 
Quantum phase flows are analyzed using the star-product technique which we believe to be the most adequate tool 
for studying the subject. 

We give definitions and recall basic features of the method of characteristics in Sect. II. 

In Sect. Ill, fundamental properties of quantum characteristics are derived. The Weyl's association rule, the 
star-product technique, and the star-functions are reviewed based on the method proposed by Stratonovich [7|. We 
show, firstly, that quantum phase fiow preserves the Moyal bracket and does not preserve the Poisson bracket in 
general. Secondly, we show that the star-product is invariant with respect to transformations of the coordinate 
system, which preserve the Moyal bracket. Thirdly, non-local laws of composition for quantum trajectories and the 
energy conservation along quantum trajectories are found in Sect. III-D. Applying the invariancc of the star-product 
with respect to change of the coordinate system (jIII.15[) and the energy conservation, we derive new equivalent 
representations of the quantum Hamilton's equations Eq. piI.19[) - pil.211) . In Sect. III-E, we derive using the star- 
product technique the semiclassical reduction of the quantum Hamilton's equations to a system of first-order ODE 
involving along with quantum trajectories their partial derivatives with respect to initial canonical variables. Finally, 



we express the phase-space Green function 
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35[ in terms of quantum characteristics and reformulate relation 



36j using the method of characteristics. 



between quantum and classical time-dependent observables 

The possibility of finding quantum trajectories and generalized Jacobi fields by solving a system of ODE gives 
practical advantages because of the existence of efficient numerical ODE integrators. It would be tempting to extend 
method of characteristics to constraint systems such as gauge theories, relativistic QMD transport models, etc. 

The skew-gradient projection method is found to be useful to formulate classical and quantum constraint dynamics 



m 



m 
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40|. In Sect. IV, we show that in classical constraint systems characteristic lines exist and the 
method of characteristics is efficient. The proof we provide does not presuppose that constraint equations can be 
solved. The phase flow is commutative with the phase flows generated by constraint functions. Characteristic lines, 
if belong to the constraint submanifold at r = 0, belong to the constraint submanifold at r > also. 

Sect. V gives description of quantum characteristics in constraint systems. Although the formalism is complete, we 
encounter unexpected difficulty to formulate simple geometric idea that quantum trajectory belongs to a constraint 
submanifold. Using tools of the analytic geometry, any idea like that requires the use of composition of functions. In 
quantum mechanics, one has to use the star-composition. This calls for a modification of usual geometric relations 
"belong", "intersect", and others. In a specific quantum- mechanical sense, the Hamiltonian and constraint functions 
can be said to remain constant along quantum trajectories, while in the usual geometric sense they obviously don't. 
The problem of visualization of relations among quantum objects in phase space is discussed in Sects. III-D and V-B. 

Conclusion summarizes results. 



II. CHARACTERISTICS IN CLASSICAL UNCONSTRAINED SYSTEMS 



The phase space of system with n degrees of freedom is parameterized by 2n canonical coordinates and momenta 
C — (9^5 ■••■'?"iPii ■■■iPn) which satisfy the Poisson bracket relations 



{e.n 



jkl 



(II.l) 



with 



-Er, 
En 

where En is the identity n x n matrix. The phase space appears as the cotangent bundle T^M" of n-dimcnsional 
configuration space R". The matrix En imparts to T^R" = R^" a skew-symmetric bilinear form. The phase space 
acquires thereby structure of symplectic space. 

In what follows, physical observables are time dependent, whereas density distributions remain constant. Such a 
picture constitutes the classical analogue of the quantum-mechanical Heisenberg picture. 

In the classical unconstrained systems, phase flow: S, ^ C ^ c(^,t), is canonical and preserves the Poisson bracket. 
The classical Hamilton's equations (II. 3|) are first-order ODE. The energy is conserved along classical trajectories 



n{0 = nc{^,T)). 

The classical Hamilton's equations (|I.3p can be rewritten as first-order PDE: 

d 



dr 



= {cn^,r),7^(c(C,T))}. 



(11.2) 

(H.3) 
(H.4) 



The phase-space trajectories can be used to solve the Liouville equation (|I.2p which is the first-order PDE. Any 
observable /(^, r) is expressed in terms of c(^, r), as indicated in Eq.(JL4|. 
Classical trajectories obey the dot-composition law: 



c'(C,n+T2) = cXc(C,ri),r2). 



(11.5) 



III. CHARACTERISTICS IN QUANTUM UNCONSTRAINED SYSTEMS 

[tI isdi 



The Stratonovich version of the Wey l's quantization and dequantization 



discussed in the next subsection and 



in more details in Refs 



Q 



14 



15|, 



M 



45|. 



A. Weyl's association rule and the star-product 

The phase-space variables ^* correspond to operators y'^ = (q^, ..., q",pi, ■■■,pn) acting in the Hilbert space, which 
obey commutation rules 

[y^J:'] =-in/'='. (III.l) 

Operators f acting in the Hilbert space admit multiplications by c- numbers and summations. The set of all operators 
constitutes a vector space. The basis of such a space can be labelled by ^*. The Weyl's basis looks like 

The association rule for a function /(^) and an operator f has the form [7|] 

J2n 



fiO = Trmm. f = J ^^f(^)'^(0- (ni.2) 



The value of /(^) can be treated as the ^-coordinate of f in the basis 23 (^), while Tr[Q5(^)f] as the scalar product of 
03 (0 and f. 

Using Eqs. pil.2p one gets an equivalent association rule 

f = /(-^^i^) exp(-^ryfey^-)|^=o. (111.3) 

The half-Fourier transform, 

f{0 = j d^Oc^^{-'-Y,eyA <q+^-\\\q-'^->, (III.4) 

provides the inverse relation. The Weyl-symmetrized functions of operators of canonical variables have representation 
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f = /(^^^^^,p|2)), (ni.5) 

where the subscripts indicate the order in which the operators act on the right. 

Given two functions /(^) = rr[05(^)f] and g(^) — rr[Q3(^)g], one can construct a third function, 

called star-product. In terms of the Poisson operator 



one has 



f{^)^9{0^n^)^M'^V)g{0- 



The star-product splits into symmetric and skew-symmetric parts, 

ih 

The skew-symmetric part is known under the name of Moyal bracket. 

The Wigner function is the Weyl's symbol of the density matrix. In the Heisenberg picture, the Wigner function re- 
mains constant W{£^^ t) ~ W{£,^ 0), whereas functions representing physical observables evolve with time in agreement 
with equation 

^f{tr) = f{^,T)^H{0. (III.7) 

This equation is the Weyl's transform of equation of motion for operators in the Heisenberg picture. It is the infinite- 
order partial differential equation (PDE). 

The series expansions of f{S,,T) over r is given by 

CXD g 

s 

where /(^) = /(^,t = 0) is the initial data function. 



Using the •-adjoint notations of Ref. 



equation (JIII.7P can be represented in the form 
|:/(^,t) = -Ad.H[f{rj,T)m. (III.9) 



Its formal solution. 



fit t) = eM~rAd,H)fit 0), (III.IO) 

is equivalent to Eq. (IIII.8[) . If the target symbol /(<^, 0) is semiclassically admissible, the evolution operator has 
asymptotic expansion [20| 

CXD 

expi-rAd^H) = ^ h^'j^^'^r). (III.ll) 

s=0 

The power series expansion in h is valid for semiclassically admissible symbols H and /. If, however, / is a rapidly 
oscillating symbol, then (jIII.lip fails and the solution of the evolution equation becomes of the WKB type whose 
exponential phase is a symplectic area (see for details Ref. i49i]). 

B. Quantum phase flow^ preserves the Moyal bracket 

Active transformations modify operators f and commute with 23 (^). Passive transformations change the basis and 
keep operators fixed. These views are equivalent. We choose the former. Consider transformations depicted by the 
diagram 



: 


u 


1 


f 


> 


y 



where il is given by Eq. (jl.ip . 

The operators of canonical variables are transformed as f —^ p = il+y'il. The coordinates ^' of new operators y' in 
the old basis 25 (^) are given by 

r ^ e = u\^, r) = Tr[^{Oii+m- (ni.12) 

Since il is the evolution operator, functions u*(^,r) can be treated as the Weyl's symbols of operators of canonical 
coordinates and momenta in the Heisenberg picture. For t — 0, we have m*(^, 0) = ^^ 

The set of operators of canonical variables is complete in the sense that any operator acting in the Hilbert space 
can be represented as a function of operators f. One can indicate it as follows: f — /(y). The Taylor expansion 
of /(y) permits the equivalent formulation of the Weyl's association rule. Transformations f ^ f — il+fil generate 
transformations of the associated phase-space functions: 

fiO^AO - /(e,T)-rr[»(Oii+fii] 



1 d^fjO) 
s\dS,'^...d^' 

s!9ei...ac* 
s!aei...ae 

/(*"(?, r)). (III.13) 






^gs#T^"'-«-) —'•«-' 



Last two lines define the star-composition. The star- function f{*u{^,T)) is a functional of w(^,r). The o-product 
is not associative in general. However, the indices is for s = l,...,2n are symmetrized, so the order in which the 
o-product is calculated is not important. 

The antisymmetrized products p^^ ...j:*^=l of even number of operators of canonical variables are c-numbers as a con- 
sequence of the commutation relations. These products are left invariant by unitary transformations: il+j:['^...j:*^''lil = 
j:['^..j:*^=l. In phase space, we get u^^^^{^,t) -k ... 7k-u'2=l(^, r) = ^['^ • ... •^*^=1 and, in particular, 

wXC,r) Aii^'(e,T) =f AC^' = ~r^. (III. 14) 

Phase-space transformations induced by il preserve the Moyal bracket and do not preserve the Poisson bracket, so 
the evolution map (, ^ £, = u{^, r), is not canonical. Using Eq. (jIII.27p . one can check e.g. that for H{(_) = {Sij^'^^^)'^ 
where 6ij is the Kronecker symbol functions m*(^, e) do not satisfy the Poisson bracket condition for canonicity to 
order 0{e^h'^). 

For real functions w*(C, t) satisfying Eas. (jIII.14p one may associate Hermitian operators y* which obey commutation 
rules for operators of canonical coordinates and momenta. As a result, functions u^{^,t) appear in the coincidence 
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with a unitary transformation relating f and p. The conservation of the Moyal bracket for a one-parameter set of 
continuous phase-space transformations is the necessary and sufficient condition for unitary character of the associated 
continuous transformations in the Hilbert space. 

C. Change of variables which leaves the star-product invariant 

Applying Eq. (jIII.13|) to product fg of two operators, we obtain function f{() *g(C)lc=*M(C,T-) associated to operator 
il'*'(f0)il and function /(*w(C,t)) -k g{-ku{£,,T)) associated to operator (il+fil)(il+gil). These operators coincide, so do 
their symbols: 

f{C)*g{C)k=*u(i,r)^f{MC,T))*g{Mtr)). (III.15) 

The star-product is calculated with respect to C, and ^ in the left- and right-hand sides, respectively. Equation (IIII.lSp 
is valid separately for symmetric and skew-symmetric parts of the star-product of the functions. 

The substantial content of Eq. piI.15[) is that one can compute the star-product in the initial coordinate system and 
change variables ^ —>■ C ^ *"(Cj ''')j or equivalently, change variables ^ -^ C — *^(Cj ''") ^^'^ compute the star-product, 
provided Eq.(lIIL14|) is fulfilled. 

The functions m*(^,t) define quantum phase flow which represents quantum deformation of classical phase flow. 




FIG. 1: Schematic presentation of the star-composition law (jIII.16[) . The solid line stands for a quantum trajectory u'(5,r) = 
ii'(*u(^, s),T ~ s) at < r < t. The dashed line is assigned to a trajectory u'(it(^, s),t — s) which we would have at s < r < i 
for the classical dot-composition law. The distance between the solid and dashed trajectories is of order of tV^. 



D. Composition law for quantum trajectories and energy conservation law 

In the usual geometric sense, quantum characteristics u(^, r) cannot be considered as trajectories along which 
physical particles move. The reason lies, in particular, in the star-composition law 

u(C, n + T-2) - u{M^, ri), ra) (III.16) 



which is distinct from m(^,ti + T2) — u{u{^,ti),T2), see Fig. [T] In classical mechanics, the composition law has the 
form of Eq. (|II.5| ). 

The energy conservation in the course of quantum evolution implies 

HiO=HiM^,T)) (111.17) 

where H(^) = Tr[^{^)^] is Hamiltonian function. H{^) is, however, not conserved along quantum trajectories in the 
usual geometric sense, so -ff(^) ^ H{u{£^,t)). In classical mechanics, the conservation law has the form pi.2[) . 

To express the idea that a point particle moves continuously along a phase-space trajectory, one has to use the 
star-composition (jIII.16|) . The dot-composition is not defined in quantum mechanics. 

Similarly, H(u{S,,t)) does not make any quantum-mechanical sense. One has to work with H{-ku{^,T)). If so, 
the only way to express quantitatively the fact of the energy conservation along a phase-space trajectory is to use 
Eq. (|i!lT71) . 

The similar problem arises in constraint systems when we want to decide if quantum trajectories belong to a 
constraint submanifold. 

The analytic geometry provides tools to formulate relations among geometric objects. Those relations which are 
expressed through composition of functions are modified. We discuss if possible to assign a geometric sense to formulas 
involving the star-composition in Sect. V-B. 

E. Reduction of quantum Hamilton's equations to a coupled system of ODE for quantum trajectories and 

generalized Jacobi fields 

Quantum Hamilton's equations can be obtained applying the Weyl's transform to evolution equations for Heisenberg 
operators of canonical coordinates and momenta 

— u'(e,r) = u\^,t)AH{0 (III.18) 

= u\tT)AH{M^,T)) (III.19) 

= CAif(C)lc=.„(5,r) (111.20) 

= {C,i^(C)}lc=*«(?,r). (111.21) 

To reach the step 2, the energy conservation piI.17[) is used. Going from (JIII.19P to pil.20p . the change of variables 
piI.15|) is performed. To achieve pil.21[) . we exploit ^' A /{£,) = {d/(C)}- The time derivative of u*(^, r) can 
be computed classically using the Poisson bracket. The substitution C = *u(^,t) leads, however, to deformation of 
classical trajectories. Equations pll.20|) and pll.21|) are the quantum analogues of Eq. (II.3p . Eq. llIII.lSP is the quantum 
analogue of Eq. pi.3p . and Eq. pil.19p is the quantum analogue of Eq. pi.4p . 

As distinct from the de Broglie-Bohm trajectories (see e.g. [43|), w*(^,t) are not related to specific states in the 
Hilbert space. 
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The functional form of quantum Hamilton's equations pil.181) is left invariant by the change of variables ^ ^ v = 
*^-(C) provided the map v^: v = w-(^), preserves the Moyal bracket. 

Equations (|III.18|) are not invariant under canonical transformations. Consider e.g. canonical map: {q,p) -^ (Q, P), 
with generating function S2{q,P) = qP + q'^ + qP^ such that p = dS2{q,P)/dq and Q — dS2{q, P)/dP. One can 
compare fog and fAg in the coordinate systems {q,p) and (Q, P). For functions f — q and g — p, one gets, respectively, 
/°.9l(,,p) = 9P / f°9\{Q,P) = qp+6h^Q/{l + 2Pf+0{h^) and /A.g|(,,p) == 1 ^ /A.g|(Q,p) = l + 2Ahy{l + 2Pf+0{h^). 
The symmetric and skew-symmetric parts of the star-product are both not invariant under canonical transformations. 
Coordinate systems in phase space if related by a canonical transformation provide non-equivalent quantum dynamics. 
This ambiguity is better known as the operator ordering problem. 

The quantum deformation of classical phase flow can be found by expanding 

oo 
s=0 

The right-hand side of Eqs. piI.18|) F^{C) = {C:H{()} is a function of (^ = •M(f,r) (i.e. functional of u(^,t)), so we 
have to expand 

oo 



using e.g. the cluster-graph method 



20 



s=0 



46|. Classical trajectories Uo(^,r) satisfy classical Hamilton's equations 



d 

— u^ = F^{uo) 



and initial conditions Uq(^,0) = ^*. Given Uq{^,t), the lowest-order quantum correction u\{^,t) can be found by 
solving first-order ordinary differential equations (ODE) 

16' "''^^'^^ °^'^'^ a^^^a^^^ 

_ _}_Tkihjk2i2 ni ji2 p3 ^ -^o("o) 

with initial conditions u\{^, 0) = 0. The functions Jq j, and Jq j,; entering Eq. (jIII.22[) is a particular case of generalized 
Jacobi fields 

Given ■u^(^,r) and J^ ^ j, (r, ^) for < ?' < s, the next corrections uI^-j^{^,t) can be found from first-order ODE 
involving generalized Jacobi fields (|III.23p with < r < s. For a harmonic oscillator, u* (^, r) = for s > 1, in which 
case quantum phase flow is both canonical and unitary. 

The generalized Jacobi fields pil.23p satisfy ODE also. The lowest order equations have the form: 

^_J^ _ dF^M .JJJ24) 



JrM...kM-)- ^ck: .1 - (ni-23) 



11 



The first of these equations describes the evolution of small perturbations along the classical trajectories. Being 
projected onto a submanifold of constant energy it becomes the Jacobi-Lcvi-Civita equation 41]. In stochastic systems, 
Jq j, grow exponentially with time. 

At any fixed level of accuracy of the semiclassical expansion, we have a coupled system of ODE for ^^.(r, ^) and 



J: 



r,ki...kt 



(r, ^) subjected to initial conditions 



«o(c,o)-r, 

<(C,0) = 0, 






(III.25) 
(III.26) 



where r > 1 and r > 1 or i > 2, respectively. The evolution problem can be solved e.g. numerically applying efficient 
ODE integrators. 

A numerical computation of the semiclassical expansion of the quantum phase ffow in the elastic scattering of 
atomic systems is performed in Ref. [29|. 

An alternative approach allowing to reduce the semiclassical quantum dynamics to a closed system of ODE is 



proposed by Bagrov with co-workers 



30 



31 



,l3a 



33|. The phase-space trajectories appearing in 



30 



M 



M 



33| are 



connected to specific quantum states like in the de Broglie - Bohm theory. 

Properties of quantum paths, localization of quantum systems, and a coherent-type representation of the quantum 
ffow are discussed in Ref. 34| . 



The series expansions of u'(^,r) and f{^,T) — /(*w(C,t)) over r are given by 

OO g 

s 

OO g 

^^(...((/(OAF(0)AiJ(e))A...fl'(e)). 

s=0 ' ' """^""^ 



(III.27) 



(III.28) 



In general, quantum phase ffow is distinct from classical phase flow. This feature holds in integrable systems also, 
as discussed in Appendix A. 

The lowest order operators 7^*-' entering Eq. ljIII.lip can be found to be 



201, 



291 



4a, M 



- ^4^(e, r)4„(e, r)4''™(C, r)/K(e, r))..,fc. 
Here, the derivatives of f{uo{£,,T)) are calculated with respect to Uq: 

5^/K(e,T)) 



(III.29) 
(III.30) 



/("oCC^t)),;!. 



Ouq ...duQ^ 



(III.31) 
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The Jacobi fields with the upper indices are defined by 

J;''=-'='(^,t) = /'-^^^■^../'-^^^^ j;^,...^,(e,r). (III.32) 

According to Ea. pil.29p . time dependence of the zero order term "/'■-^^t) f (S^) is determined by time dependence 
of the classical phase-space trajectory and form of the function /(^). The similar conclusion holds for 7^^''(''')/(C)- 
Ea. pil.30p tells that time dependence enters through the classical trajectory, the first quantum correction to the 
classical trajectory, and the classical Jacobi fields. The problem of convergence of a formal power series expansion is 
always a difficult subject. The convergence rate of the time series depends obviously on the system. Generally, the 
series expansion in r has a finite convergence radius. However, in the itegrable systems one has a truncated series 
expansion (cf. Eq.(III.27) and Eq.(A.3)). 

F. Green function in phase space and quantum characteristics 

Using orthogonality condition 



and Eq. (|III.13|l . we express Green function for the Weyl's symbols 21 



, |35| in terms of the quantum characteristics: 



= i2nhr6''\Mtr)-0 

= (27rft)"(52"(^-*-a(C,-T)). (III.33) 

A compact operator relation between the classical and quantum time-dependent observables is established in Rcf. 



36l |. Solutions of the quantum and classical Liouville equations, /(^, t) and /c(^,''"), with initial conditions /(C,0) = 
/c(C, 0) are related through the product DD^^ where Dc is the classical Green function 

Z),(C, C, r) = (27r?i)"<52"(c(e, r) - C). (III.34) 

In terms of the characteristics, we obtain 

/(?, r) = /c(c(*u(e, r), -r), r). (III.35) 

It is assumed that classical and quantum hamiltonian functions coincide i.e. %{£,) = H{^). 
Given the Green function is known, the quantum trajectories can be found from equation 

-^(C,r)=y (^m^,C-r). (III.36) 

For H = 1 — jr^e where e is an infinitesimal parameter, the associated transformations of canonical variables and 
phase-space functions are given by S£,^ = ^* A e-ff(0 = {C'l^-^IO} ^^^ ^fiO — fiO ^ ^^(0- The transformations 
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of canonical variables are canonical to order 0(e) only. The infinitesimal transformations of symbols of operators are 
not canonical. Any function H{S,) can be used to generate classical phase flow or quantum phase flow, according as 
the dot-product or the star-product stands for multiplication operation in the set of phase-space functions. 

The analogue between unitary and canonical transformations is illustrated by Dirac [17| in terms of the generating 
function S{q',q) defined by exp(j-S{q' ,q)) =< q'\ii\q >■ The evolution map {q,p) —^ {q',p'), is canonical for p = 
—dS{q', q)/dq and p' — dS{q' , q)/dq' . The parallelism of the transformations is manifest, but trajectories are complex. 
The generating function defined by the phase of < (7'|il|g > yields real trajectories. It is not clear, however, if time- 
dependent symbols of operators are entirely determined by such trajectories. 

The Weyl's symbols of operators of canonical variables u'(f,r) are the genuine characteristics in the sense that 
they allow by equation /(^, r) — /(t^w(^, r), 0) the entire determination of the evolution of observables. The quantum 
dynamics is totally contained in m*(^,t), whereas the deformation of symbols of the operators calculated at *M(^,r) 
has a kinematical meaning. 

IV. CHARACTERISTICS IN CLASSICAL CONSTRAINT SYSTEMS 

We give first description of second-class constraints systems and of the skew-gradient projection formalism. The 
details are found elsewhere [IJ, [ij, Isvl, ba ba, |40 1 . 

A. Classical constraint systems in phase space 

Second-class constraints Ga{0 = with a ~ 1, ..., 2777, and m < n have the Poisson bracket relations which form a 
non-degenerate 2777 x 2777 matrix 

det{g,(O,et(O}^0. (IV.l) 

If this would not be the case, it could mean that gauge degrees of freedom appear in the system. After imposing 
gauge-fixing conditions, we could arrive at the inequality pV.l[) . Alternatively, breaking the condition pV.ip could 
mean that constraint functions are dependent. After removing redundant constraints, we arrive at the inequality 
(lIVTl) . 

Constraint functions are equivalent if they describe the same constraint submanifold. Within this class one can 
make transformations without changing dynamics. 

For arbitrary point f of the constraint submanifold F* = {^ : Ga{0 = 0}j there is a neighbourhood where one may 
find equivalent constraint functions in terms of which the Poisson bracket relations look like 

{GaiO.Gbim^^ab (IV.2) 
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where 



Here, Em is the identity m x m matrix, lablb, 



E„ 

~Em 

The matrix I''^ 



(IV.3) 



-Tnh is used to hft indices a, 6, 



up. 



The basis (|IV.2p always exists locahy, i.e., in a finite nei 



This is on the fine with the Darboux's theorem (see e.g. 



jhbourhood of any point of the constraint submanifold. 



4l|). AU symplectic spaces are locaUy indistinguishable. 



idGb(e 



IdGaCn 




FIG. 2; Schematic presentation of skew-gradient projection onto constraint submanifold along commuting phase flows generated 
by constraint functions. 



B. Skew-gradient projection formalism 



The concept of the skew-gradient projection S,s{S,) of canonical variables ^ onto a constraint submanifold plays 
inrportant role in the Moyal quantization of constraint systems. Geometrically, the skew-gradient projection acts 
along phase flows IdQ°-{^) generated by constraint functions. These flows are commutative in virtue of Eqs. (|IV.2p : 
Using Eqs. pV.2[) and the Jacobi identity, one gets {^°, {Q^, /}} — {Q'', {Q"^, /}} for any function /, so the intersection 
point with F* is unique. 

To construct the skew-gradient projections, we start from equations 



{UO,Ga{0}=0 



(IV.4) 



which say that point ^s{£,) S F* is left invariant by phase flows generated by Ga{0- Using the symplectic basis (jIV.2[) 
for the constraints and expanding 

1 



UO = ^ + X'^Ga + ^X'^'GaGb 



(IV.5) 



in the power series of Qa^ one gets 



CXD ^ 

xGaiGa2---Ga^- 



k=0 



(IV.6) 
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Similar projection can be made for function f{£,): 

oo _. 

xGa^Ga^-Ga,- (IV.7) 



fc! 

fe=0 



One has 

fsiO = fiUO)- (IV.8) 

The projected functions are in involution with the constraint functions: 

{fsiO,GaiO}^0- (IV.9) 

Consequently, /s(^) does not vary along IdQa{S), since 

{/(C),5(c)}^^(/d5(c)r. 

The skew-gradient projection is depicted schematically in Fig. [2l 

C. Evolution and skew-gradient projection 

In the classical second-class constraints systems, one has to start from constructing Ws(C) from 7i(^). The evolution 
equation for phase-space functions can be converted then to the classical Liouville equation: 

^/(C,r) = {/(?, r),H.(C)} (IV.IO) 
Similarly, the canonical variables obey the classical Hamilton's equations: 

^c*(e,r) = {c^(e,r),H.(0} (IV.U) 
with initial conditions 

c'(C,0)=r. (IV. 12) 
Equation 

{ea(O,W.(?)}-0 (IV.13) 

tells that Qa{£,) remain constant along dlS^jT): 

Qa{0=Qa{c{S„T)). (IV.14) 

Equations pV.14p show that trajectories do not leave level sets {^ : Qa{C) ~ constant} and therefore do not leave 
the constraint submanifold F* = {^ : Qa{0 = 0}. 
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Given Ti.s{£,) is constructed, it becomes possible to extend standard theorems of the Hamiltonian formahsm to 
second-class constraints systems without modifications. The novel element is the interplay between the evolution and 
the skew-gradient projection. 

Let the coordinate system {ry*} is obtained from the coordinate system {^*} by the canonical transformation 
^ ^ 77 = c(^,r). 



Ck.^ phase 



flow 



c(e,T) 




FIG. 3: Classical phase flow c{(,,t) is commutative with classical projection £,3{^) onto constraint submanifold F*. 



Eq. (jIV.7|) may be applied for c^{(^,t). Using Eq. (jIV.14p . we replace the arguments of the constraint functions 
to ^{^jT) and replace everywhere d{^,T) with ?/*, as long as the Poisson brackets are invariant and the constraint 
functions are scalars. We arrive at 



= 6(c(^,t)). 



(IV. 15) 



The first line is a consequence of Eq. (jIV.8[) . The evolution is commutative with the skew-gradient projection. Equation 
PV.ISP is illustrated on Fig. H 

The Liouville equation can be solved provided phase-space trajectories c(^, r) are known. In general, 



Applying projection pV.7p . one gets 



/(e,r)=/(c(C,r),0). 

/.(e,T) = /(c(6(C),t),0) 
= /(c,(^,r),0). 



(IV. 16) 



(IV. 17) 



The first line follows from Eg. pV.Sp . Equation pV.17p shows how to use characteristics in order to solve evolution 
equations in the classical second-class constraint systems. 

The evolution depends on choice of the constraint functions up to a canonical transformation. Suppose we found 
two sets of the constraint functions Ga{C) ^nd GaiO describing the same constraint submanifold. Each set can be 
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transformed to the standard basis (|IV.2|) . Such bases are related by canonical transformations, so one can find a 
canonical map: ^ — > v = w-(^), such that Gai^) = Ga{v-{^)). The inverse transform is u ^ ^ = t;-|_(u). The 
skew-gradient projections ^^(0 ^^^id Vs{v) are related by: 

Vs{v-{0)^v^{UO)- (IV.18) 

The skew-gradient projection depends on choice of the constraint functions up to a canonical transformation. The 
same is true for projected Hamiltonian functions: 

nsiO = Kiv) (IV. 19) 

where ^.'{v) ~ 7i(u_|_(u)). Two sets of the constraint functions GaiO a-nd Ga{0 l^ad to the canonically equivalent 
Hamiltonian phase flows. 

V. CHARACTERISTICS IN QUANTUM CONSTRAINT SYSTEMS 

The constraint systems represent high interest since all fundamental interactions in the elementary particle physics 
are based on the principles of gauge invariance. Gauge fixing turns gauge-invariant systems into constraint systems. 

In the classical mechanics, the constraint systems can be treated as a limiting case A ^ oo of systems in a potential 
V\{q) which rapidly increases when the coordinates q go away from the constraint submanifold. In the limit of 
A — > oo, Vx{q) = if g belongs to the constraint submanifold and V\{q) = +oo when q does not belong to the 
constraint submanifold. The classical systems obtained by imposing the constraints and by the limiting procedure 



have equivalent dynamic properties |4l| . 







In the quantum mechanics, this is not the case. The limiting procedure applied to a particular system of Ref. |42| to 
model holonomic constraints, results to the quantum dynamics which depends on the way the constraint submanifold 
is embedded into the configuration space. From other hand, the quantization of constraint holonomic systems leads to 
the conclusion that the dynamics is determined by the induced metric tensor only 



401 



43| . The limiting procedure and 



imposing the constraints are not equivalent schemes of the quantization. In what follows, we discuss the constraint 
dynamics as it appears in the gauge theories. 

The Groenewold-Moyal constraint dynamics has many features in common with the classical constraint dynamics. 
The projection formalism developed for constraint systems allows, from other hand, to treat unconstrained and 
constraint systems essentially on the same footing. 

A. Skew-gradient projection in quantum mechanics 

We recall that classical Hamiltonian function 7i(^) and constraint functions Ga{0 a-i'^ distinct in general from their 
quantum analogues H{^) and Ga{£.)- These dissimilarities are connected to ambiguities in quantization of classical 
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systems. It is required only 

lim H{0 = niO, lim G,(e) - GaiO- 
In what follows 

r* = {e : GaiO = 0}. (V.l) 

The quantum constraint functions Ga{^) satisfy 

Ga{0/\Gl,{0=Iab. (V.2) 

The quantum-mechanical version of the skew-gradient projections is defined with the use of the Moyal bracket 

6(OAGa(O=0. (V.3) 

The projected canonical variables have the form 

fc=0 

oGa,oGa,...oGa,. (V.4) 

The quantum analogue of Eq. (jIV.7P is 

oo - 
oGa,oGa,...oGa,. (V.5) 

The function ftiO obeys equation 

/t(C)AGa(e)=0. (V.6) 

The evolution equation which is the analogue of Eq. (JIII.7P takes the form 

§if{0 = f{0/\Ht{0 (V.7) 



where Ht{£,) is the Hamiltonian function projected onto the constraint submanifold as prescribed by Eq. (|V.5p . 

Any function projected quantum-mcchanically onto the constraint submanifold can be represented in the form [l5| 

MO = ^Kt(C))- (V.8) 

In the space of projected functions, the set of projected canonical variables £,t{£,) is therefore complete. 
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FIG. 4: Quantum projection ^tiO defined by Eq.(4.21). The submanifold F* = {Ct(?) : (. G T^R"} does not coincide with 
the constraint submanifold F* = {^ : Ga{£,) ~ 0}. The variance is of order ~ h . The constraint submanifold F* can be 
parameterized by classical projection F* = {^s(0 : C G T»R"} constructed with the use of the quantum constraint functions 

B. Coordinate star-transformations do not keep geometric relations among quantum objects 



The evolution equation in the quantum constraint systems has the form of Eq. (|V.7p which is essentially the same 
as in the quantum unconstrained systems. Replacing H{^) by Ht{^), one can work further with solutions u{^,t) of 
quantum Hamilton's equations (|III.18[) . It is not required for points ^ to belong to the constraint submanifold, so 
phase-space trajectories u{^, r) occupy the whole phase space. 

The quantum phase flow preserves the constraint functions in the following sense: 



GaiO^GaiM^T)). 



(V.9) 



The alternative equation Ga{^) — Ga{u{^,T)) which would carry the conventional geometric meaning uses pre- 
conditionally the dot-composition law which is not allowed quantum-mechanically. It is obviously violated, so in the 
usual sense u{^, t) ^ T* for r > even if u{^, r = 0) = ^ £ T* (see Fig. [5]). 

Any attempt to decide if w(^, r) G T* involves the dot-composition e.g. 



u{^,t) er* ^ VaGa(u(C,T))=0. 



(V.IO) 



Statements involving the dot-composition are, however, forbidden. 

Surprisingly, expressive means of the star-product formalism are not enough to formulate the simple geometric idea 
that a trajectory belongs to a submanifold. 

We wish to find statements admissible quantum-mechanically and from other hand which would support relations 
of belonging and intersection inherent for geometric objects. 
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It is tempting to interpret Eas. (|V.9p as the evidence that quantum trajectories u(^,t) do not leave, in a specific 
quantum- mechanical sense, level sets of constraint functions {f : Ga{S.) — constant}. 

Such a statement has the invariant meaning with respect to unitary transformations: Suppose the map u+: v —>■ (, = 
v+{v), corresponds to a unitary transformation in the Hilbert space. The inverse unitary transformation generates 
the inverse map u_: £, ~* v — W-(C), such that u_(*w_|_('L')) — v and, by virtue of Eg. piI.lSp . v+(*v_(^)) — ^. In the 
coordinate system {i^*}, the constraint functions become 

G'M = Ga{*v+iv)). (V.U) 

Equation (jIII.15|) allows to change the variables f — > -kv+{v) in Eq. (IV.9p to give 

G'Jv) = G'J*u'iv,T)) (V.12) 

where 

u'{v,t) ^v_{*u{*v+{v),t)) (V.13) 

represents the quantum phase flow in the coordinate system {w*}. Equations (jV.9|) and (|V.12p are therefore equivalent. 
They show that "do not leave" represents a predicate invariant under unitary transformations. 

The non-local character of relations between the quantum phase flows is displayed in Eq. (|V.13|) explicitly. One can 
conclude that quantum trajectories do not transform under unitary transformations as geometric objects. 

u'{v,t). 




FIG. 5: Constraint submanifolds F* and F*' (solid lines) and quantum trajectories u{(^,t) and u'{v,t) (dashed lines) in unitary 
equivalent coordinate systems {^'j and {v'}, respectively. As shown, u{^,t) crosses F* twice, whereas its image u'{v,t) crosses 
F" once. Any counting of the intersections rests on an implicit use of the dot-composition, an operation which is forbidden 
quantum-mechanically. The property of the statements u(^,r) G F* and u'{v,t) £ F*' be true or false depends on unitary 
transformations. From the viewpoints of Eqs. (|V.9|) and (|V.12|) . u{^, r) and u'{v, r) belong to the level sets of Ga(^) and G'a{v), 
respectively. However, from condition Ga(*u(^,r)) = it does not follow that G'a(*u' {v , t)) — and vice versa. Geometric 
relations among quantum objects, which use the dot-composition, do not have objective meaning. 

The coordinate transformation 7j_ : ^ — > w = «_ (^) does not superpose T* and 

r' = {v : G'M = 0}. (V.14) 
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Assuming ^ £ T*, we obtain G'^{v_{£,)) ^ G'^{*v_{£,)) = Ga{0 = and therefore w_(^) ^ T*' in general. The 
constraint submanifold does not transform under unitary transformations as a geometric object also. 

We see that points of F* transform differently from F*. They are "not attached to F*". In new coordinate system, 
F* represents a set of new points. To put it precisely, 

^eF*^w = w_(0 eF*' = w_(F*). (V.15) 

Unitary transformations affect the visualization of trajectories and submanifolds. The relation "do not leave" 
supports, however, some features inherent to the usual geometric relations "belong" and "intersect". One can show 
e.g. that if quantum trajectories do not leave the level sets of Ga{£,) and each level set of Ga{£,) is a subset of one of 
the level sets of Fa{£,) then quantum trajectories do not leave the level sets of Fa{^). 

One cannot assign to quantum trajectories definite values of energy and constraint functions. In the coordinate 
system {^'} one has E^ ~ H{iru{£^,T)), whereas in the coordinate system {v^} one has E^j — H'{-ku'{v,T)) where 
H'iy) is defined by Eq. (jA.2[) . The constants E(^ and E^ do not depend on time. However, E^ ^ E^ in general even 
if trajectories are related by a unitary transformation. The same conclusion holds for constraint functions, as shown 
on Fig. O 

Finally, the syntax of the star-product formalism is not rich enough to express the simple geometric idea that 
trajectory belongs to a submanifold. 

The star-product geometry admits the statement that quantum trajectories do not leave level sets of the constraint 
functions. The validity of this statement is not affected by unitary transformations and has the objective meaning. 
The quantum- mechanical relation "do not leave" is the remnant of usual relations of belonging and intersection 
inherent to geometric objects. It cannot be completely visualized, however. 

C. Evolution and skew-gradient projection 

The classical phase flow commutes with the classical skew-gradient projection, as discussed in Sect. IV. We want 
to clarify if such a property holds for quantum systems. 

Given the quantum trajectories u(^, r) are constructed, the evolution of arbitrary function can be found with the 
help of Eq. pil.i3p and its projection can be computed using Ea. (|V.5[) . 

The quantum projection applied to arbitrary function cannot be expressed in terms of the same function of the 
projected arguments Eq. (jV.8[) . basically because the classical relation {fg)s = fs9s turns to the quantum inequality 
if * 9)t 7^ ft * 9t- In terms of a function ip{^) defined for /(^) = /(^,0) in Eq. (jV.8[) . the quantum analogue for 
Eqs. (|IV.17P reads 

M^,r)^^{*Ut{^,T)). (V.16) 

The construction of (/?(^) from /(^) is a complicated task, so practical advantages of this equation are not seen 
immediately. 
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Equation (jV.16p accomplishes solution of the evolution problem for observable /(C, t) in terms of quantum charac- 
teristics. 

It remains to prove 






(V.17) 



The first line is a consequence of the fact that the constraint functions Ga{£,) are Moyal commutative with the 
projected Hamiltonian function Ht{^). To arrive at the second line, it is sufficient to use Eq. (|V.9|) to replace arguments 
of the constraint functions entering the skew-gradient projection. 

The quantum phase flow commutes with the quantum projection, as illustrated on Fig. [6l 

The composition law (|III.16p for quantum phase flow holds for the constraint systems. It holds for projected 
quantum trajectories also: 



as a consequence of Eas. (jV.17p . 



Uti^.Tl +T2) = Uf(*Ut(^,Tl),r2), 



(V.18) 



u(e,T) 




FIG. 6; Quantum phase flow is commutative with quantum projection operation: ti(*Ct(0 >''')) ~ ^t{*u{^,T)). The phase-space 
trajectory itf(^,r) does not belong to the submanifold F* = {^t(0 • C G T,R"} except for r = 0, so the white planes on Figs. 
|4] and [6] are distinct. 



VI. CONCLUSIONS 



The method of characteristics for solving evolution equations in classical and quantum, unconstrained and con- 
strained systems has been discussed. The analysis rests on the Groenewold-Moyal star-product technique. 

The classical method of characteristics applies to first-order PDE and consists in finding characteristics which are 
solutions of first-order ODE. For the classical Liouville equation, the corresponding first-order ODE are the Hamilton's 
equations and the characteristics of interest are the classical phase-space trajectories. 
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TABLE I: Solutions of evolution equations for functions (second column) and projected functions (third column) of classical 
systems (first row) and quantum systems (second row) in terms of characteristics. c(^, r) are solutions of classical Hamilton's 
equations with hamiltonian function H(^) (second column) and projected hamiltonian function T~l3{^) (third column). Cs{(,,t) 
are classical projections of c{^,t). u(^,r) are solutions of quantum Hamilton's equations with Hamiltonian function H{^) 
(second column) and projected Hamiltonian function Ht{(,) (third column). ut{^, r) are quantum projections of u(^, r). (p{^, 0) 
is defined in terms of ft{£., 0) by Ea. (|V.16p . Classical and quantum projections are defined by Eas. pV.7[) and (|V.5|l . respectively. 



Systems: 


unconstrained 


constrained 


classical 


/WC,r),o) 


ficsii,r),0) 


quantum 


/(*M(5,r),0) 


^(*Mt(^,r),0) 



The quantum Liouvillc equation is the infinite-order PDE. Nevertheless, it can be solved in terms of quantum 
characteristics which are solutions of the quantum Hamilton's equations. These equations represent infinite-order 
PDE also. 

Using the star-product formalism, we showed that to any fixed order in the Planck's constant, quantum character- 
istics can be constructed by solving a closed system of ODE for quantum trajectories and generalized Jacobi fields. 
The quantum evolution becomes local in an extended phase space with new dimensions ascribed to generalized Jacobi 
fields. This statement holds for constraint systems also. 

One-parameter continuous groups of unitary transformations in quantum theory represent the quantum deformation 
of one-parameter continuous groups of canonical transformations in classical theory. Quantum phase flow, induced 
by the evolution in the Hilbert space, does not satisfy the condition for canonicity and preserves the Moyal bracket 
rather than the Poisson bracket. The knowledge of quantum phase flow allows to reconstruct quantum dynamics. 

The results reported in this work are valid for semiclassically admissible functions, i.e. for functions regular in h 
a.t h ^ Q. Physical observables are normally associated with classical devices and expressed as classical functions of 
classical variables. The quantum evolution turns, however, the set of classical functions into the set of semiclassically 
admissible functions. 

The use of the skew-gradient projection formalism allows to treat unconstrained and constraint systems essentially 
on the same footing. We showed that the skew-gradient projections of solutions of the quantum Hamilton's equations 
onto the constraint submanifold comprise the complete information on c^uantum dynamics of constraint systems. 

The formalism we developed applies in particular to the dynamics of gauge-invariant systems which become second 
class upon gauge fixing. The quantum dynamics of charged particles in external gauge fields on flat and curved 
manifolds is discussed within the star-product formalism in a gauge-invariant manner in Refs. 



53|, 



54| 



The evolution equations for semiclassically admissible functions admit solutions in terms of characteristics in all 
physical systems, as summarized in Table HI 

The analytic geometry uses the dot-product and rests on classical ideas how to arrange composition of functions. 
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It is well known that all theorems of geometry can be reformulated using tools of the analytic geometry. 

Given the dot-product is replaced with the star-product, we arrive at the star- product geometry with well defined 
coordinate systems, transformations of the coordinates and equations for functions of the coordinates. However, 
objects of the star-product geometry, defined algebraically, can hardly be visualized: 

We found that quantum trajectories and constraint submanifolds do not transform as geometric objects. The 
statement "quantum trajectory belongs to a constraint submanifold" can be changed to the opposite by a unitary 
transformation. The star-composition law (jIII.16|) shows also that the quantum evolution cannot be treated literally 
as moving along a quantum trajectory. 

We attempted to find statements whose validity cannot be reverted by transformations of the coordinate system 
and which, from other hand, express relations similar to "belong", "intersect", etc. A weak but consistent geometric 
meaning can be attributed to the statement "quantum trajectories do not leave level sets of constraint functions". 

The dot-product composition of linear functions coincides with the star-product composition of linear functions, so 
under linear transformations straight lines and hyperplanes turn to straight lines and hyperplanes. Relations of the 
linear algebra, imbedded into the star-product geometry, preserve the consistent geometric meaning. 

Finally, this work extended the method of characteristics to quantum unconstrained and constraint systems. From 
the point of view of applications, it is motivated by the fact of using classical phase-space trajectories in transport 
models and by the appearance of constraints in relativistic versions of QMD transport models. The method of quantum 
characteristics represents the promising tool for solving numerically many-body potential scattering problems. 
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APPENDIX A: QUANTUM PHASE FLOW IN INTEGRABLE SYSTEMS 

A completely integrable classical system admits a canonical transformation which makes the Hamiltonian function 
depending on half of the canonical variables only (see e.g. [41j). Such variables if exist can be taken to be canonical 
momenta which usually referred to as actions. The canonically conjugate coordinates are referred to as angles. 

In quantum mechanics, we search for a unitary transformation (or a half-unitary transformation [l8|) allowing to 
express the Hamiltonian as an operator function of operators of canonical momenta. If such a transformation exists, 
the Hamiltonian commutes with the operators of canonical momenta, so that the canonical momenta are integrals of 
motion, whereas the operators of canonical coordinates depend linearly on time. 

The quantum integrable systems admit an equivalent treatment in the framework of the Groenewold-Moyal dy- 
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namics 55|. Given the Hamiltonian function H{^) is known, one has to search for a map w+: 

v^C = v+{v), (A.l) 

preserving the Moyal bracket, for which the system admits a Hamiltonian function 

H'{v)^H{-kv+{v)) (A.2) 

depending on actions 1;"+"'^, . . . ,t'^", i.e., canonical momenta only. We restrict the discussion by unitary transforma- 
tions (jA.ip . leaving aside more involved cases described in Refs. [18J, ll9|. 

Let m*(^,t) and o^^u.t) be solutions of Eq. ljIII.lSp with Hamiltonian functions H{£) and H'{v), respectively. In 
the coordinate system {f*}, the series expansion (JHI.27P is truncated at s = 1. The quantum Hamilton's equations 
give a 'motion by inertia': 

a'{v,T)=v' + {v\H'iv)}T. (A.3) 

The actions (i = n + 1, ..., 2n) remain constant, whereas the angles (i = 1, ..., n) evolve linearly with time. Equations 
(|A.3p can be derived as the Weyl's transform of the equations of motion for the Heisenberg operators of the canon- 
ical coordinates and momenta obtained by a unitary transformation from the initial set of operators the canonical 
coordinates and momenta. 

The Poisson bracket {w*, H'{v)} depends for any i on the actions only, so one has 

a'(w, r) o a"'(u,T) = a^{v,T)a^ {v,t), 
a'{v,T) Aa^{v,T) = {a*(w, r), a-' (w, r)}. 

The map a: 

V ^ V ~ a{v,T), (A. 4) 

showing the evolution in the coordinate system {u*} is both canonical and unitary, as the left-hand side of Eq. (|A.3p 
is a first-order polynomial with respect to the angles. 

The actions u"+^, . . . , u^" Poisson and Moyal commute with H'{v). Composite functions wL (•u(*w_|_(i;), t)), where 
V- is the inverse unitary map: 

^^v = v.{0, (A.5) 

such that u!_(*w+(w)) — w*, obey Egs. pn.lSp and proper initial conditions and coincide with a''{v,T). It can be 
expressed as follows: 

u\(„t)^v\{M*v^{0,t)). (A.6) 

The functions v± are defined using the star-product and depend on h accordingly. 
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We thus conclude that quantum phase flow is distinct from classical phase flow for integrable systems also. For 
a one-dimensional system H = ^p^ + V{q), which is a classical integrable sysem for any potential V{q), the first 
quantum correction to the phase-space trajectories appears to order 0{h^T^). 

In general case, Ea. (|A.6|) shows the connection between quantum phase flows u*(^,t) and a^{v,T) in two unitary 
equivalent coordinate systems {^'} and {v^}. 
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